Are there significant scalar resonances in B DW°K-K° decays? 
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Abstract 

We study the indication of large different branching ratios between B~ — > d(*'°K~K° and 
Bd — » K~ K° observed by Belle. Interestingly, the same situation is not found in the decays 
B —* D(*> K~ K*° . If there exist no intermediate resonances in the decays B~ -> D^°K~K , 
a puzzle will be arisen. We find that the color-suppressed processes B~ -* D^°a (U50) with 

(1450) — ► K~K° could be one of the candidates to enhance the BRs of B — >■ D^°K~K°. Our 
conjecture can be examined by the Dalitz plot technique and the analysis of angular dependence 
on K~K° state at B factories. 
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Since CP violation was discovered in .fT-meson system in 1964 [lj, the same thing has 
been also realized by Belle j^J and Babar P with high accuracy in the B system. Although 
the new observation doesn't improve our recognition on CP violation, it stimulates the 
development. Besides the origin of CP violation, B factories also provide the chance to 
understand or search the uncertain states, such as the scalar bosons / (400 — 1200), / (980), 
ao (980) and glueball etc. Unlike pseudoscalar mesons, the scalar bosons with wide widths 
are difficult to measure directly via two-body B meson decays. Inevitably, the study of three 
body decays will become important for extracting the quasi-two-body resonant states from 
the Dalitz plot and the invariant mass distribution. By the analyses of Dalitz plot technique, 
we can recognize whether there exist unusual structures in the phase space. 

Recently, Belle has observed the BRs of the decays B -> D^R~R*° to be BR(B~ -> 
DW°K-K*°) = 7.5 ± 1.3 ± 1.1 (15.3 ± 3.1 ± 2.) x 10~ 4 and BR(B° -> D^ + R~R*°) = 
8.8 ± 1.1 ± 1.5 (12.9 ± 2.2 ± 2.5) x 10" 4 , and the decays B -> D^R~R° to be BR(B~ -> 
DW°K-K°) = 5.5±1.4±0.8 (5.2±2.7±1.2) x 10" 4 and BR(B° -> D^ + R~R°) = 1.6±0.8± 
0.3 (2.0 ± 1.5 ± 0.4) x 10" 4 [4]. According to Belle's analyses, we know that R~R*° system 
has the state J p = 1 + ; and also, it is pointed out that the B — > D^*' R~ R*° decays in the low 
R~R*° invariant mass region can be fitted well by quasi-two-body decays B — > D^a^ (1260) 
with a^(1260) -> R-R*°. In R~R° system of B~ -> D°R-R° decays, by adopting the 
fitting distribution dN/dcos9KK oc (Rl cos 2 6kk + (1 — Rl) sin 2 9kk) Belle observed the 
value Rl = 0.97±0.08, i.e. the R~R° state prefers being J p = 1~. Interestingly, by looking 
at the data shown above, we find that the BRs of charged B decaying to R~R° final state are 
much larger than those of neutral B decays. In terms of central values, the ratios, defined by 
RM = BR(B- -> D^R-R°)/BR(B° -> D^+R'R ), could be estimated roughly to be 
3.5 (2.6); however, there is no such kind of implication on the R~R*° final state. Therefore, if 
the data display the correct behavior, there must be something happened in B — ► D^R~R° 
decays. Before discussing the origin of the differences, we need to understand the mechanism 
to produce R~R*° and R~R° systems in charged and neutral B decays. 

Since B — > D^K~K^*'° decays are dictated by the b — > cud transition, we describe the 
effective Hamiltonian as 

H cS = ^K[Ci(//)0i + C 2 (//)(9 2 ] (1) 
with Oi = daiipCpba and = daUaCpbpi where q a % = Qa 7^(1 ~ IhjQpi a (P) are the color 




FIG. 1: Topologies for B — * D^K K^*'° decays, (a) color-allowed with q = u and d, and (b) 
color-suppressed. 

indices, V u = V* d Vch are the products of the CKM matrix elements, Ci^tjit) are the Wilson 
coefficients (WCs) p. As usual, the effective WCs of a 2 = Ci + C 2 /N c and a\ = C 2 + Ci/N c , 
with N c = 3 being color number, are more useful. According to the interactions of Eq. 

we classify the topologies for B — > D^K~K^° decays to be color-allowed (CA) and 
color-suppressed (CS) processes, illustrated by Fig. [TJ and also use the decay amplitudes 
Afi [K~ K^°) to describe their decays, where Q = C(N) denote the charged (neutral) B 
decays. From the figure, we see that both CA and CS topologies contribute to charged 
B decays but neutral B decays are only governed by CA. We note that because we don't 
consider the problem of direct CP asymmetry, the annihilation effects, which can contribute 
to Bd — > D^*^ K~ decays and are usually smaller than the CA contributions, are 
neglected. Although Fig. ^ corresponds to the picture of three-body-decay, after removing 
the ss pair, the figures are related to quasi-two-body decays if du could form a possible 
bound state. 

Now, we have to examine that if there exists only J p = 1~ state in K~K° system, 
the differences between BR(B~ -> D^ K-K ) and BR(B d -> D^+R-R ) are insignif- 
icant. To be more understanding the problem, it is useful to start from the discussion 
on B -> D^R-R° processes. As mentioned before, charged B decays are governed 
by CA and CS while neutral B decays are only from CA; therefore, the amplitudes for 
B~ and B° decays can be written as A C (R-R°) oc aiA^(A^A"°) + a 2 M% (R~R°) and 
A N (R~R°) oc d\M.a {R~ R°), where ■M^n t \ express the hadronic transition matrix ele- 
ments for CA (CS). For simplicity, we only consider factorizable parts. As a conjecture, if 
there exist significant differences in BRs of B~ and B° decays, the source should be from 
a 2 M.fr(R~R°). In the following, we use two examples to show that it is impossible to 
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only increase the BRs of B -> D^K~K° without enhancing the BRs of B -> D^K~K*°. 
Firstly, if the decays Z? — > D^*'K~K° are only governed by quasi-two-body decays, say 
B — > D(*>px with px being arbitrary vector meson, the decay amplitudes can be described 
by A C (K~K°) ~ a i/ w F B - DW (0) + a 2 / flW F B ^(0) and ^(#-#°) ~ o^F^^O), 
where / D (.) and / px are the decay constants of and px mesons, respectively and 
jpB-tD *) {px) d erio te the B — > D^*\px) form factors. However, it is known that \a2/a\\ ~ 0.26 
and f Px F B ^ DM (0) ~ f DM F B ^P x (0) 6]. If the resonance is J p = 1~ state, the dif- 
ferences in BR between charged and neutral B decays will be slight. The same argu- 
ment can be applied to B — > D^K~K*° decays for J p = 1 + resonant state. Hence, 
the consequences are consistent with the observations of B — > D^*' K~ K*° decays. Sec- 
ondly, we consider nonresonant mechanism on B — ► D^K~K°. Since the situation cor- 
responds to a three-body phase space, we use the K~K° invariant mass, expressed by 
u, as the variable to describe the behavior of the decay. Hence, the decay amplitudes 
could be written as A c (K~ K°) ~ ai F^ KK (u 2 )F B ^ DW (u 2 ) + a 2 f DW F B ^ KK (u 2 ) and 
A N (K~K°) ~ ai F°^ KK (uj 2 )F b ^ dM (u 2 ), and the p°^ KK (F B ^ KK ) are the associated form 
factors. It is known that at asymptotic region, in terms of perturbative QCD (PQCD) the 
behavior of F°~^ KK has the power law 1/uj 2 (1iiu 2 /A'q CD )~ 1 pj. If we take the behavior of 
p*B->D * £ k e — uj 2 /M x ) n , Mx is the pole mass and n = 1 or 2, we clearly see that 
the decay amplitudes are suppressed at large u. Furthermore, according to the analysis of 
Ref. [a], the dominant region actually is around u 2 ~ Am^. Therefore, if a 2 F B ^ KK (Am^) 
is the source of the differences in BR between charged and neutral B decays, the simi- 
lar effects a 2 F B ~* KK * (Am B ) will also affect the decays B — > D^K~K*°. Hence, unless 
F B ~^ * (Auib) is much less than F B ^ KK (Arris), otherwise, the CS effects should have the 
similar contributions to B — > D^*'K~K*°. It is known that with the concept of two-meson 
wave functions the calculations of F B ^ KK (Arris) and F B ~* KK * (Arris) can be associated 
with the wave functions of KK and KK* systems. At uj 2 ~ ArriB region, two- meson wave 
functions could be described by two individual meson wave functions 101. Refer to the 
derived K and K* wave functions the value of F B ^ KK (Arris) / F B ^ KK * (Arris) should 
be close to F B ~* K (0) / F B ^ K * (0) ~ 0.8. As a result, it is hard to imagine that the CS effects 
on pure three-body decays can make large differences in B — > D^K^K processes but not 
in B — ► D^K'K* processes. 

Inspired by Belle's observation, we think the significant differences in the BRs of B — > 
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D^K^K should not come from K~K° with J p = 1~ state but with J p = + state. This 
could be easily understood as follows: as discussed before, only CS effects, a<iM^ \K~ K°), 
could make a discrepancy in BRs of charged and neutral B decays. Therefore, the mechanism 
which dominantly contributes to CS topology will be our candidate. If K~K° has the J p = 
+ state, due to F°^ KK oc < KK, + \u n f ll d\0 >, in terms of equation of motion, we get that 
po-*KK - g p r0 p 0r tional to (m^— m u ). That is, the contribution of scalar state to CA topology 
is negligible. On the contrary, there is no such suppression to the CS topology. Now, the 
oroblem becomes how to produce K~K° to be a + state. By searching particle data group 



12| . we find that the preference could be the scalar bosons ao(980) and ao(1450) because 



both are isovector states and have sizable decay rates for a — > KK. Hence, to realize our 
thought, we propose that the decays B~ -> D^°a (980, 1450) with (980, 1450) -> K~K° 
can satisfy the required criterion to enhance the BRs of B~ — ► D^*'° K~ K° . 

Since the quark contents of scalar mesons below or near 1 GeV are still obscure in the 

n 

literature for avoiding the difficulty in estimation, we only make the explicit calculations 
on a (1450) which have definite composed structure of qq. If regarding a (980) consists of qq 
state, the same estimation could be also applied Q|. In theory, it is known that the serious 
problem on two-body nonleptonic B decays comes from the calculations of hadronic matrix 
elements. Since the involving processes are governed by CS topologies, like the well known 
decays B -> J/#.K"W and B d -> D«° 7T° in which nonfactorizable effects play an important 
role, we adopt the PQCD approach which is described by the convolution of hard amplitude 



and wave functions 



can deal with the factorizable and nonfactorizable contributions 



and can avoid the suffering end-point singularities self-consistently [16]. 

In the calculations of hadronic effects, the problem is how to determine the wave functions 
of and ao mesons. For £)M meson wave functions, we could model them to fit with the 
measured BRs of Bd — > D^°tc° decays 13], in which B meson wave function is chosen to 
be coincide with the observed BRs of B — ► tttt decays while 7r meson ones are adopted from 
the derivation of QCD sum rules ^l|. As to ao scalar meson, the spin structures of a are 
required to satisfy (0|-u'y /i (i|a ( 7 ,p 3 ) = [(m d - m u )/m ao }fp% and (0|gg|a ) = m a J, where m ao 
and / are the mass and decay constant of ao- In order to satisfy these local current matrix 
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elements, the distribution amplitude for clq is adopted as 



dxe ixP - z {V\mjq{z)i\ao) 



V2N. 



(2) 



By the charge parity mvanance and neglecting the effects oflight current quark mass m <d) , 

— x) [181 ] . and their normalizations are 



-$ O0 (l-x) and $P Q (x) 



■ a ' 



we obtain $ ao (x) 

J" dx$ ao (x) = and J Q dx^^ Q (x) = f/2\/2N c . Although vector D*° meson carries the spin 
degrees of freedom, in the B~ — > D*°a^ decay only longitudinal polarization is involved. 
The results should be similar to D Therefore, we only present the representative 

formulas for B~ — > D^a^ decay. Hence, the decay amplitude for B~ — > D°a^ is read as 



A = V,, 



DO 



where J :e { y M e ^) is the factorized (non-factorized) emission hard amplitudes. According to 
Eq. (j2J), the typical hard functions are expressed as 



D°a 



M D o aQ 



7] J dx\dx 3 / bidbib 3 db 3 ^B(xi,bi) 
o Jo 

(Ci + ( 2 x 3 )$ ao (x 3 ) + r«(l - 2x 3 )$ p aQ (x 3 ) 
S e (tl)+r a (2^ ao (x 3 )-r a ^ ao (x 3 )) 

rl poo 

2rj / d[x] / b 1 db 1 b 3 db 3 ^ B (x 1 ,b 1 )^ D (x 2 ) 
Jo Jo 

{ [ - ((1 - r 2 a )x 2 + (1 - r 2 D )x 3 ^ ao (x 3 ) 
+r a x 3 ^ ao (x 3 )] #(tj) + [((Ci - r a 2 )(l - x 2 ) 
+r 2 a x 3 - r 2 D x 2 )§ ao {x 3 ) - r a x 3 $ p aQ (x 3 ) 



(3) 



(4) 



with V = 8-kC f M 2 , r a{D) =m ao{D) /M B , d = l-r 2 a -r 2 D , ( 2 = Ci~r 2 D , = a.(t*) a 2 (4) 

fiW„(f e ) &,({*},{&}) and = a s (f d )(C7 2 (^)/iV c ) Su(t d ) B+D+ao h d ({x},{b}). tg, Su 

and h et d denote the hard scales of B decays, Sudakov factors and hard functions which are 
arisen from the propagators of gluon and internal valence quark, respectively. Their explicit 
expressions can be found in Ref. jl^ ]. 



For numerical estimations, the B, and ao meson wave functions are simply chosen 



as 



$ B {x,b) = N B x 2 (l - x) 2 exp 
3 



[-51 


'xm B s 


f co 2 B bh 




\ UJ B > 


1 2 J 



*SoO«0 



2^2^ 



2V2iVc 



/ DW x(l-x)[l + 0.7(1 -2x)], 
Qx(l - x)Cl ,2 {\ - 2x) 
3(1 -2x) 2 ~\ 



(5) 



where N B is determined by the normalization of dx& B (x, 0) = f B /(2^2N c ) and C^ 2 (y) is 
the Gegenbauer polynomial. The values of theoretical inputs are set as: uj b = 0.4, f B = 0.19, 
f Di *) = 0.20 (0.22), / = 0.20, m B = 5.28 and m DM = 1.87(2.01) GeV. By the taken values 
and using Eq. (jHJ) with excluding WC of a 2 , we immediately obtain the B — > a (1450) 
form factor at large recoil to be 0.44. The result is quite close to the value 0.46 which 
is estimated by the light-cone sum rules (20]. Hence, the magnitudes of the considered 
hard functions are given in Table UJ Consequently, the BRs of B~ — > /}(*)°a~(1450) are 

TABLE I: Hard functions (in units of 10~ 2 ) for B~ -> D^°a Q (U50) decays with uj b = 0.4, 
f B = 0.19, f DM = 0.2 (0.22) and / = 0.20 GeV. 





M D°a 


fD*^ B *o aQ 




-1.42 


-2.22 + i0.97 


-1.56 


-2.39 + »1.06 



obtained to be 8.21 (9.34) x 10" 4 . Since the predictions of PQCD on the BRs of B -> L>°7r° 
and B -> J/^K^ and the helicity amplitudes of B -> J/^iT 1 are consistent 
with the current experimental data, with the same approach, our results should be reliable. 
Furthermore, the BR products of Br(B- -> D^a^ (1450)) x Br{a^ (1450) -> JsTfT ) « 
1.81 (2.05) x 10~ 4 with Br(aQ (1450) -> if-fT ) ~ 0.22 [12]. Clearly, the contributions of 
quasi-two-body decays to B~ — > D^°K~ K° modes are close to the pure three-body decays 



5 D^ + K~K° 



22 



In summary, we have investigated that when a proper scalar meson is considered, the 
BRs of Br -> D^°K~K° will deviate from those of £ d -> D^+K~K°. Although we only 
concentrate on ao(1450), the same discussion is also applicable to ao(980). Since our purpose 
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is just to display the importance of scalar boson on the decays B — > D^*'°K~K°, we don't 
consider the theoretical uncertainties at this stage. And also, we neglect discussing the 
interfering effects of resonance and non-resonance. It is worthful to mention that by using 
powerful Dalitz plot technique, many scalar mesons in charm decays have been observed 
by the experiments at CLEO H, E791 H, FOCUS Q, and Babar 2(J- Expectably, 
the significant evidences of scalar productions should be also observed at B factories. In 
addition, since the scalar meson doesn't carry spin degrees of freedom, there is no specific 
direction for K~K° production so that we should see a different angular distribution of 
K~K°, such as the coefficient associated with the term sm 2 9 kk in the fitting distribution 
mentioned early will be enhanced. Hence, with more data accumulated, our conjecture can 
be examined by the Dalitz plot and the analysis of angular dependence on K~K° state. 
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